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Abstract 

We propose an alternative mechanism for the gating of biological membrane channels in response 
to membrane tension that involves a change in the slope of the membrane near the channel. Under 
biological membrane tensions we show that the energy difference between the closed (tilted) and 
open (untilted) states can far exceed k^T and is comparable to what is available under simple di- 
lational gating. Recent experiments demonstrate that membrane leaflet asymmetries (spontaneous 
curvature) can strong effect the gating of some channels. Such a phenomenon would be more easy 
to explain under gating-by-tilt, given its novel intrinsic sensitivity to such asymmetry. 



The correct biophysical function of mechanically sensitive membrane channels, such as the 
widely studied MscL[l-4], is vital in maintaining the viabihty of hving cells. These channels 
are typically made up of a 5-8 transmembrane proteins that form a barrel-like assembly[5]. 
Fluid flows through a central pore in the open (active) state but is either restricted, or 
suppressed entirely, in the closed (inactive) state. As well as maintaining osmotic balance 
mechanosenstive channels play important sensing roles in touch, hearing, turgor control 
in plant cells etc [6-11]. The fact that MscL functions in reconstituted membranes [2] is 
good evidence for a gating mechanism that depends only on membrane tension, rather than 
cytoskeletal effects [12] or signalling cascades. We propose a new mechanism for the gating of 
biological membrane channels in response to elevated membrane tensions, which we will refer 
to throughout as gating- by-tilt, see Fig 1. The elevated membrane interfacial tension that 
acts to open the channel can be generated in several ways including via an osmotic imbalance 
between the cell interior and exterior or by changes in the cells morphology during adhesion, 
fihpodia formation etc[13]. 

The gating mechanism of such mechanosensitive channels have been probed by numerous 
patch clamp experiments (see e.g. [7] and references therein) and have been studied fairly 
extensively by molecular dynamics (MD) simulations, see e.g. [14]. The results of these 
studies have been interpreted as being consistent with gating that operates primarily via a 
process resembling simple dilational opening of the channel, see Fig 1. However, both the 
experimental techniques and MD simulations have limitations that make it difficult to be 
conclusive concerning this mechanism. Current computational limits on the MD simulations 
allow only for only very short simulation times, typically of the order of a few nanoseconds, 
while it is known that gating takes milliseconds to occur after a step change in the patch 
clamp pressure that controls membrane tension[7]. It is common practice in MD simulations 
to artificially apply thermodynamically large forces to the membrane (channel) in order to 
induce it to open on the available timescales. It is therefore not clear that these simulations 
are able to differentiate between certain models for how these channels gate. Indeed it 
is one goal of the present work to propose a new paradigm for consideration in future 
simulations. Patch clamp experiments essentially measure the tension at which the channel 
opens and the size of its open pore by way of conductivity measurements[7, 15, 16]. It isn't 
clear how these experiments can be used to differentiate between different gating models. 
Interestingly it is also known that the originally proposed pentameric structure of MscL [5] 
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may be somewhat different from its in vivo structure [7]. Thus it is possibile that some 
membrane tilt may be induced by the orientation of the walls of mechanosensitive channels 
in their closed configuration. Such a feature must be present if gating-by-tilt is to play a role. 
During the opening of the channel the membrane tension does work to change the channel's 
conformation. The conformational free energy of the channel itself (and alone) changes by 
an amount AGchannei = Gopen — ^closed which is positive, indicating that the closed state 
has a lower conformational free energy. In what follows we will address only whether the 
channels are open or closed at equilibrium under constant membrane tension. The existing 
paradigm for tension-gated channels involves a simple dilational transition from the closed 
to the open state with an associated increase in total channel area AA [7]. This is the 
increase in the total effective "footprint" of the channel within the membrane. The change 
in free energy is proportional to this and to the membrane tension 7 



and is the only energy that is available to overcome the change in conformational free energy 
of the channel AGchannci- The proposed gating-by-tilt is driven by changes in the slope of 
the membrane where it meets the channel. The simplest version of this involves a constant 
membrane slope 9 in the closed state which relaxes {9 0) after opening, see Fig 1. A 
central pore could be opened in association with such a transition. In practice a hybrid 
mechanism may operate in which there is both some change in tilt and some dilation. Other 
workers have attempted to study how changes in the lateral pressure profile of the channel 
might also affect gating[17-19] although, in what foffows, we wiif restrict our attention to 
effects that depend directly on the membrane tension. 

We argue that it is possible for the action of the membrane tension 7 to efficiently gate 
the channel provided the work done by the membrane AGmemb is large enough 



where AGmemb = AGdiiatc or AGtiit, or a combination of the two. The criteria of Eq (2) 
equate to a channel that can be substantially closed below tension 7 and substantially open 
above it. It is reasonable to assume that nature has been able to evolve a channel with a 
configurational energy change that is roughly optimal for gating at the desired tension, in 
which case one expects AGmemb ~ AGchannei- If AGjaemb/kgT is larger than unity then a 



AGdilate = 



(1) 



AGinemb(7) ^ AGchannel k, kgT 



(2) 
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substantial channel configurational energy can be overcome by the action of the membrane 
tension and hence the channel can remain substantially closed at low tensions and substan- 
tially open at high tensions (above the gating threshold). Alternatively if AGmemh/kgT is 
smaller than unity the energy available from the membrane tension is inadequate to over- 
come any activation barrier that is greater than the thermal energy scale and so the channel 
is either predominantly closed, if AGchannei > ^s^^; or opens and closes rather randomly, 
if AGchannei < k^T. Ftom Eq (2), and the arguments above, it is possible for a channel 
that opens by gating-by-tilt to efficiently gate the channel at tension 7 when the parame- 
ter controlling this efficiency AGtiit{'y)/kgT > 1. We now proceed to calculate the energy 
AGtiit. 

The Hamiltonian for membrane that is asymptotically fiat, with normal parallel to the 
2;-axis, but which has a small normal deviation from flatness, of magnitude ^(r) due, e.g. to 
the presence of a membrane channel, is given by [20] 

H^lj dh[K{v\r+^{vur] (3) 

where k is the membrane rigidity, r is the radial position with r = |r| (see Fig 1) and V 
is the two dimensional (plane polar) version of the operator. The total energy associated 
with distortion of the membrane can be established by a variational approach on Eq (3) and 
depends on the boundary conditions for the membrane. Up to unimportant global rotations 
or translations of the entire frame the displacement of membranes which are asymptotically 
flat at infinity is found to be 

u = aKQ(kr) + /3\ogkr (4) 

for r > a where a and (3 are yet undetermined constants, k = yj^j 1^ is an inverse length 
characteristic of the membrane and Kq is the usual modified Bessel function of the first 
kind of order zero. It can be shown that solutions of the form of w ~ log kr correspond 
only to channels which exert a finite integrated normal force on the membrane [24] , e.g. by 
anchoring onto the cytoskeleton or an external substrate. Thus for channels which exert 
no overall normal force we have u — aK^ikr) where the constant a is fixed by a boundary 
condition corresponding to the angular tilt at the periphery of the channel Vii(r — a) — 
—ak Ki{ka) — —9. Hence 

u = eKo{kr)/(k Ki(ka)) forr > a (5) 
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where the resulting normal membrane deviation (a few nm or less) falls to zero over a 
distance ~ from the edge of the channel.The free energy difference AGtiit between a 
channel in the closed state and one in the open state {9 0, say) follows by substitution of 
Eq (5) into Eq (3) integrated over r > a. 
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AGtiit = TTukae'-p—^ (6) 
Ki[ka) 

A quantitative comparison with AGdiiate is shown in Fig 2. 

For steady state physiological tensions in the range 7 = 10~^-10~^N/m one has 0.03 < 
ka < 0.1 and the system is quantitatively inside the regime ka va. which AG^iit can be 
shown to have the following analytic approximation 



AGtiit ~ 7r7a'^2 



2 

log Euler Gamma 

ka 



forA;a<l (7) 



Where Euler Gamma has a constant value of about 0.6. At the highest gating tensions 
ka can approach unity. Interestingly Eq 1 and Eq 7 appear to be very similar, both being 
the product of an area change and the tension 7 except for the factor of log ^ appearing 
in Eq (6). In order to test whether gating- by-tilt is indeed a plausible mechanism that 
may contribute to AGmemb we use data collected for various Mechanosenstive channels [15] , 
from which estimates for the the area change AA, gating tension and conformational energy 
change AG have been obtained, assuming dilational gating. This gives us enough informa- 
tion to calculate the corresponding value for AGtiit from Eq 6, assuming instead a purely 
gat ing-by- tilt scheme that occurs at the same tension. Since this energy will depend on the 
angle 9 a sensible approach seems to be to ask instead what value of 9 would be required 
to entirely account for the free energy change AG (previously estimated assuming dilational 
gating) and also what angle would be required to account for (say) 10% of it. We will 
denote these angle ^100 and ^10 respectively. If a channel opens, via gating- by-tilt, to expose 
a pore of area /S.A = irb^ then we impose the boundary condition that the membrane slope 
is 9 in the closed state at radial distance r = a = d + b/2 from the central (z) symmetry 
axis of the pore, where d is the thickness of the "walls" of the channel, see Fig 1. This 
is a somewhat arbitrary, although reasonable, choice for a which is merely the radial dis- 
tance to the exterior of channel in the closed state. For the six channels for which data 
is available [15] we find (taking d = 2nm) that 39° > ^100 > 16°, while 12° > ^10 > 5° is 
correspondingly smaller. We find that the smallest values of these angles both correspond 



to the channel known as MscMJ which gates at the lowest value of AGtiit although not the 
lowest membrane tension. This may therefore represent the best candidate for future MD or 
experimental studies on channels that gate by tilt. The above angles are (all) rather small 
and thus we have demonstrated that a small change in the tilt of the membrane during 
gating can provide a significant contribution to the free energy change AGraemh- 

It has also been observed that the composition of lipids in the membrane can have a 
significant effect on the gating of MscL channels [17]. It is possible to construct arguments 
that can explain a dependence on membrane composition in terms of mismatches between 
the membrane (thickness) and the channel geometry, or the pressure distribution exerted 
by the membrane interior on the edge of the channel. However the dramatic effect recently 
observed on addition of "conical" lysophophatidylcholine (LPC) lipids to one leaflet of the 
membrane [17] are generically difficult to explain within such models. In these experiments 
it was observed that the channels will open, even under small applied tensions, only if the 
conical lipids are localized in one leaflet, giving rise to an asymmetry that tends to make 
the membrane prefer to bend into a convex shape, away from the LPC rich face, rather than 
the opposite concave one. It is possible to analyze this difference within a model that has 
the LPC homogeneously distributed across the membrane or, experimentally, the membrane 
patch that is clamped. In this case the Hamiltonian for the membrane includes the so-called 
spontaneous curvature Cq, which is a function of the difference in LPC concentrations across 
the inner and outer leaflets of the membrane. If the LPC is locahsed in the upper leaflet 
Co < 0. When the spontaneous curvature is small <^ k'^ one obtains [20] 



The corresponding minimum energy solution Eq (5) for the displacement is unchanged under 
these conditions. However the total free energy of the membrane can be shown to be 
AGmemb = AGtilt + ^Gco with AGtiit as before (6) and 



This free energy difference between the open and closed state can be positive, favoring 
opening, or negative, helping to maintain a closed channel if ^ > or ^ < respectively. 
The sign of 9 depends on whether the the channel is oriented either "up" or "down". In 
experiments on reconstituted membranes it is very likely that the channels will be in both 
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orientations. Eq (9) suggests that for moderate spontaneous curvatures Cg — l/(10a) the 
"up" channels (say) are driven to open by an energy difference A.Gc„ ~ T/c^T, even in the 
absence of any tension. This gives a Boltzmann weight of about 1000, with a corresponding 
enhancement in the number of open channels. This could help to explain why some channels 
isolated in a patch clamped membrane are observed to remain open under these conditions, 
even under small tensions (pressure differences) [17]. Indeed these authors independently 
suggest that, "The asymmetry of the lateral pressure profile between the two leaflets of the 
bilayer is what actually initiates the sequence of mechanical transduction steps that leads to 
the open state" . We argue that such effects are highly suggestive of a role of the gating-by- 
tilt mechanism proposed here. In conclusion the energy difference available between between 
the closed (tilted) and open (untilted) channel conformations is found to be comparable to 
that for dilational opening, as estimated from experimental measurements, and could be even 
greater for artificially engineered membrane systems. Thus we propose that the gating-by-tilt 
mechanism should be considered a viable candidate for channel gating, particularly in view 
of the lack of precise knowledge concerning the channel geometries in the open and closed 
states. There seems no obvious a priori reason why gating-by-tilt should not be commonly 
employed in nature, short of some generic difficulty in bio- engineering a channel architecture 
that gates-by-tilt. Furthermore, under gating by tilt, the tilting walls of the channel can 
be attached by a hinge to a fairly rigid frame, from which they swing open. This suggests 
that gating-by-tilt may have further inherent design advantages over dilational gating in 
which the entire channel must move (dilate). Since preparing this work we have learned 
of an independent study [22] that investigates the gating of mechanosensitive channels via 
effective line tensions. These would also appear within our model for large ka ^ 1, although 
we make no assumptions about the value of this parameter, which rarely seems to greatly 
exceed unity. Furthermore, these authors seem to assume that any channel tilt remains fixed 
throughout channel opening, in which case the tilt uniformly favors the closed (undilated) 
state. This is fundamentally different to the present work in which we propose tilt variation 
as the gating mechanism. 
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FIG. 1: Idealised sketch of a cylindrically symmetrical membrane channel showing two schemes for 
tension mediated gating. The upper two images show the closed and open states of the channel 
under the dilational gating model and the lower two images the same for gating-by-tilt. Under 
both schemes the membrane tension does work by increasing the combined projected area of the 
membrane and channel. 
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FIG. 2: The work available from the coupling of the membrane tension to a conformational change 
of a membrane channel AGmemb is plotted against the logarithm of the membrane tension 7. The 
dashed lines corresponding to the energy change AG^iiate available from dilational opening (1) 
and the solid lines the energy AGtiit available from gating-by-tilt (6). Three curves are shown 
in each case corresponding to different values of the channel size. For dilational gating (dashed 
curves) the radius of the central pore in the open configuration is, from upper curve to lower curve, 
6 = 1.5, 1, 0.5 nm, while for gating-by-tilt the radial distance of the exterior of the closed channel 
is, from upper curve to lower curve, a = 4, 3, 2nm. These represent biologically reasonable sizes 
[3, 15]. In this illustrative example we have taken the change in tilt to be ^ = 7r/6 = 30° and 
the membrane rigidity k = 20A;^T. The viability of the gating process relies on the value of the 
parameter AGmemb/^s^ being larger than unity, as discussed in the text. Steady state biological 
membrane tension are in the 10-^N/m to lO'^N/ m range [21] while Msc channels can gate at 
tensions up to > 10~^N/m [15]. 
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